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1. Introduction 

The dissipative harmonic oscillator has for long attracted considerable interest as a 
prototype of open quantum system. The early work of Magalinskii pQ and Ullersma [2] 
focusing mainly on an Ohmic environment and the weak coupling limit has been later 
extended in many different aspects [SI H El El 13 El El CHI CD 112 EE d as, for instance, 
to strong coupling [TT] . to non-Markovian noise [HI or to nonlinear coupling j^j. 
Much activity was boosted by the work of Caldeira and Leggett ^H] on the dissipative 
mechanics of a macroscopic quantum variable. Quantum decoherence has also been 
addressed in Refs. IT7| ITH] . A comprehensive review can be found in the textbook 
by Weiss |Hj|. 

In the above context relatively little attention has been paid to effects arising from 
the coupling to different system variables. The system variable which couples to the 
heat bath is most often assumed to be the position q. In the sequel we refer to this model 
as the g-oscillator. The choice of q as the coupling variable was favoured in Ref. j!5j . 
There it was argued that a complex dissipative environment can be modelled by a bath of 
harmonic oscillators, with the couplings parameters chosen to yield a Langevin equation 
for the q variable in the semiclassical limit. It was shown later [20] (see also Ref. [2*T] ) 
that, in a superconducting Josepshon juntion, the oscillator bath model can be derived 
microscopically from the coupling of the phase variable to the quasiparticle bath, with 
the phase playing the role of position. At that time the possibility of having a second 
bath coupled to the momentum p variable was ruled out. However, in a superconducting 
weak link, the electromagnetic field couples to the relative number variable, which in 
the above scheme plays the role of momentum [2H E2|- If the coupling to the phase is 
neglected, the semiclassical behavior of the particle number is governed by an Abraham- 
Lorentz equation [2H ES| • The coupling to the momentum has been also been considered 
by Leggett [24J and by Cuccoli et al. [25J. A systematic study of the combined role of 
the electromagnetic and quasiparticle fields is presented in Ref. [21]. Thus Josephson 
junctions provide a scenario where the distinction between the different baths and their 
different coupling mechanisms becomes imperious. 

Our study has also been motivated by the work of Castro Neto, Novais and 
coworkers |2Sj- To investigate the problem of a spin | impurity in an antiferromagnetic 
environment, they studied a generalized spin-boson problem using renormalization 
group techniques. The localized spin couples with two components to the environment. 
They found a crossover from incoherent to coherent transitions in the spin correlation 
functions in certain parameter ranges. There the term quantum frustration was coined, 
which refers to the reduction of the effective interaction with the environment due 
to the non-commutativity of the coupled operators. In this article, we argue that a 
harmonic oscillator coupled through position and momentum to two independent baths 
is a suitable representation of a large spin impurity in a ferromagnetic system. We find 
a weaker cancellation of the effects of the two baths. Since that situation corresponds 
to a large (quasiclassical) spin, we have referred to this weaker form of frustration as 
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quasiclassical frustration [27] . 

Here we wish to investigate systematically the scenario where a particle couples 
through its position and momentum to two baths, which are independent but have the 
same temperature. As a paradigmatic system we choose a harmonic oscillator. This 
choice is motivated not only by mathematical convenience. A more fundamental reason 
is that the symmetrical roles played by position and momentum makes the harmonic 
oscillator the natural ground to study the difference and the interplay (when both are 
present) between position and momentum coupling. 

An explanatory comment is in order. As is well known, a coupling to q can be 
converted into a coupling to p through a suitable canonical transformation. However, 
the two representations bear an important difference. Only in one of them is it 
possible to represent the system-bath interaction as the coupling of the quantum 
variable to a system of otherwise independent oscillators. A characteristic example 
is provided by a particle of charge e in an electromagnetic field, where a "velocity- 
coupling model" [2H] seems to apply. Minimal coupling (p — > p — eA/c) generates not 
only an interaction term p ■ A but also a diamagnetic term cx A 2 which can be rightly 
interpreted as an interaction between the effective oscillators. A unitary transformation 
U = exp(ieq • A/c), with [q il pj]_ = iSij, acting on p — (e/c)A, removes the coupling 
p • A. This happens at the expense of generating a coupling oc q • E between the position 
and the electric field E. Importantly, in this new representation no quadratic field term 
is left, i.e. the charge couples to a set of independent photons. Thus, in the precise 
language which we propose here, a charged particle couples to the electromagnetic field 
through its position q. 

In Sec. ED we study the model of an oscillator coupled to two different baths. The 
general formalism is sketched in Sec. 12.11 In Sec. 12.21 we investigate the case where the 
two independent Ohmic baths couple linearly to position and momentum. A related 
feature is the possibility of inducing a crossover from overdamped to underdamped 
oscillations by increasing the damping coefficient. Surprisingly, we obtain underdamped 
equilibrium oscillations for arbitrarily strong coupling. Sec. 12.31 is devoted to the 
properties of spectral functions under rather general combinations of double-bath 
systems. Time dependent phenomena and implications for quantum decoherence are 
discussed in Sec. Eland Sec. EJfor the case of a double Ohmic bath. In Sec. Elwe consider 
the situation where an oscillator couples to a single bath with the most general form 
of linear coupling to both position and momentum. We find that the most prominent 
feature of the double-bath model, namely, the existence of underdamped oscillations for 
arbitrarily strong coupling, disappears in the case of single-bath dissipation. 
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2. Two independent baths 

The general form of a Hamiltonian describing an oscillator coupled to two independent 
baths is [2T] 

H = ^{p + 5p) 2 + ^{q + 5q) 2 + ^u k a\ k a qk + ^u k a ] pk a pk , (1) 

k k 

where = i and all operators are dimensionless and h — 1. The form of the 

Hamiltonian highlights the symmetry between q and p. The notion of mass is avoided 
by m = 1/ujp. For the fluctuating pieces we assume that they are linear in the bath 
variables and independent of p and q, 

5q = i^2^k (a ] pk - a pk ^ 

k 

Sp = iJ2 X k Kfc - a Qk) > ( 2 ) 



k 



After two unitary transformations U p = exp(ip5q) and U q = exp(iq5p) one arrives at 
the Hamiltonian 



\ 2 

H = —q + 2_^u k a qk + —q + — p + ^ 



k 



u k - 2 



k 



l^k 

a p k H p 

u k 



2 



(3) 



where the short-hand notation a^a = \a\ 2 is used. The model is a Caldeira-Leggett- 
type-of Hamiltonian ^Hj- It describes a harmonic oscillator with momentum p and 
position q, each variable being coupled to a different oscillator bath. The frequency of 

1 /2 

the central oscillator is uq = (u p u q ) 1 . The baths are described by the spectral densities 

J g (u) = 2^2\\ k \ 2 5(u - u k ) , J p (u) = 2^ \fi k \ 2 5(u - u k ) ■ (4) 

Although one could have started directly with Eq. (j3J), we prefer to present the 
Hamiltonian Eq. as the starting point in order to provide a natural justification for 
the renormalization terms ^2 k \nk\ 2uJ k ~ p 2 and J2 k (Afcl 2 ^ 1 q 2 in Eq. @ which otherwise 
would have to be introduced ad hoc. 

For low-lying excitations, and thus for low temperatures, the Hamiltonian (jH)) 
becomes equivalent to that of a large spin s in a ferromagnetic environment. This 
can be seen as follows: The Hamiltonian for a spin in a large magnetic field along the z 
direction is 

H' = -^ B (S z B z + SJB x + S y 5B y ) + B 2 J (5) 

where the fluctuating terms model the low-lying bosonic (magnon) excitations of the 
ferromagnet at the site of the large spin. The fluctuations in the z direction are 
neglected since they are quadratic in S x and S y with S x , S y <C S z . In other words 
S z is approximately a constant of motion [2H|. For s — > oo the first term in (jSJ) is a 
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harmonic oscillator S z = h~\a\ 2 — h~s. The action of Si on the eigenstates of this harmonic 
oscillator is 

S z \n) = {n — s)\n) 

S+\ n ) = h\/1s — n\Jn + l\n + 1) 

« hy/2s^n + l\n + 1} = hy/2sa^ \n) 

3-\n) = hy/2s -n + l^/n\n - 1) 

~ h~\/2s^/n\n — 1) = hy/2sa \n) , 

and [S'_,S' + ] = 2/?, 2 s. This amounts to keeping only the leading order term in the 
Holstein-Primakoff transform of S + , S- jHO]- The rescaled Hamiltonian H' /2s becomes 
formally identical to Eq. (J3J) in the limit s — > oo, with ivo = ^bB z /2s and the fluctuating 
pieces scaling as 5B X = y/s~Yl ^k{a q k + a qk ) and 5B y = \/sY^^k{ a pk + a Pk )- Such large 
spins have been observed in magnetic particles pi] . 

In the following we assume a power-law behavior at u = for both spectral densities 
J q {uj)= 2^gLU a " /(uJpl~ 7t) and J p (lu) = 2<y p u ap / (cu^' 1 it) . We introduce the "phononic" 
frequency cj p h and the dimensionless coupling constants 7„, with n = q,p. Moreover 
we introduce a cutoff frequency fl n for both baths. Unless stated otherwise, we assume 

Wph = fig = fip = fi. 

2.1. General Results 

Elimination of the bath variables yields the Heisenberg equations of motion for q and p 
q(t) = u p p(t) + j ds K p (t - s)p(s) + F p (t) 

-pit) = u q q{t) + J ds K q (t - s)q{s) + F q (t). (6) 
The response kernel is defined as 

f°° J (uj) 

K n (t) = / — cos(ut)duj , n = q,p (7) 

Jo u 

and the force operator F q {t) = ^ ^k^qk exp(— io>fct)+H.c., with F p (t) defined accordingly 
In Fourier space, Eq. (0) reads 

J q (u)-uj q q + iup = F q (uj) (8) 



-xujq 



J p {u) -u p p = F p (u) , (9) 



where J n (cu) is the symmetrized Riemann transform 

M = T y/P- ^-»sgn(a;)|/(H). (10) 

J uj' (or - uj 2 ) 2 

The oscillation modes are given by the zeros of the function 

X _1 (^) = u% - uj 2 - uJ q J v {u) - UpJ q (u) + J q (u) J p (uj) (11) 
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where x( u ) is the generalized susceptibility. 

For the two baths at the same temperature the symmetrized correlation functions 
for the position Cqq\t) = M[q(t), q(0)] + ) and the momentum C pp \t) = ^([p(t),p(0)] + ) 
are obtained from xi 00 ) as follows 

C%\t) = - / lxM| 2 cosM coth(/^/2) 
71 Jo 




The corresponding expression for C pp (t) is obtained by interchanging p <-» q in 
Eq. ()12j) . The antisymmetrized correlation functions for position and momentum 
cL\t) = ([n(*),n(°)]-> (with 

n = QjP) are obtained from Eq. ()12|) in the standard 
way [HI, i.e. by substituting sin(o;t) for cos(uit) coth(/3a;/2). 

We wish to emphasize that, within the dissipation model studied here, the equations 
of motion (JHJ) are symmetric in q and p for arbitrary coupling strength. Thus the 
situation is different from that of a g-oscillator when described within the rotating- 
wave approximation, which is known to yield a model symmetric in q and p but which 
holds only in the weak coupling limit [38] . 



2.2. Ohmic coupling 

First we focus on the important special case that both spectral densities are Ohmic, a q 
= a p = 1 . Then the real parts of J q (u) and J p (u) vanish to lowest order in f2 _1 . The 
susceptibility has poles at the roots of the quadratic polynomial [5] 

ul - i{u q ^ p + uj p -f q )u - (1 + 7g7 P V 2 , (13) 

which are 



— (-i K ±VT^?) ee -ir^iC, (14) 



UJ 

u± = 

(1 + Iqlp) 

where 

ft EE W + (15) 

2o; (l + 7?7p) 

The solutions of Eq. (|13j) are either purely imaginary or a pair of complex conjugates 
depending on whether k is greater or smaller than 1. Thus, 

k < 1 Criterion A (16) 

is the commonly accepted criterion to distinguish between underdamped and 
overdamped oscillations. The underdamped region lies in a stripe of width A = 
An (1 + T] 4 ) 1 ^ 2 , with i] ee {ujq/ujp) 1 / 2 limited by the graphs of the functions 



(17) 
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In Fig.[T]the stripes of under damped oscillations, marked in the (7^, j p ) plane are plotted 
for three different values of rj. For large n = muo, corresponding to a large mass of the 
central oscillator, the stripe of underdamped oscillations becomes smaller. In the limit 
77 — > 00, the introduction of an infinesimal coupling to the momentum can induce a 
transition from overdamped to underdamped oscillations. However the range of values 
of 7 P allowing for underdamped oscillations also becomes increasingly small as oc rf 1 . 

The behavior of the system as a function of 7 P for fixed 7 g is also interesting. We 
set 7] = 1. The inverse damping time r _1 is a monotonously increasing function of 7 P 
for 7 g < 2 and a monotonously decreasing function for r y q > 2. Therefore an additional 
bath coupling to p has opposite effects on the damping time r depending on whether the 
original g-oscillator starts in the underdamped (7^ < 2) or in the overdamped (7^ > 2) 
regime. In the first case, the additional bath always reduces the damping time, leading 
to infinitely strong damping in the limit 7 P — > 00. However, in the regime 7, > 2 we can 
drive the system from the overdamped into the underdamped regime by increasing the 
coupling strength j p . What is more, we can take 7 P = 7 9 = 7, which corresponds to a 
completely symmetric Hamiltonian in p and q. At this point the system is always in the 
underdamped regime and for large 7 the oscillator frequency is close to its maximum. 
We find for the inverse damping time 



Therefore we are led to the paradoxical situation that the higher the friction coefficient 7 
the larger r becomes. In particular, for 7 — >• 00 one gets uj t — > 00. This time dilatation 
is in itself a remarkable effect, since it contrasts with the behaviour of the ^-oscillator, 
for which r — > as 7 9 — > 00 [see Eq. 114))] . 

In the general, non-symmetric case, we note that for a given finite 7 9 the damping 
time r cannot be made to acquire arbitrary values by tuning 7 P . Rather, it is bounded 
between 2^ q /u)Q and 2/7^0, which is higher depending on 7 g ^ 1. 



1 



u>o7 




T 



1 + 7 2 ' 




Figure 1. Left: Stripes of underdamped oscillations in the (7 g ,7 P ) plane for three 

1 /2 

different values of the parameter 77 = (u! q /iv p ) — mu>o, according to criterion A [Eq. 
I16| . Right: Regions of underdamped oscillations with criterions B [Eq. Q19[l] and C 
[Eq. (EH], for 17= 1. 
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However, from these striking observations one should not conclude that in the limit 
of infinite 7 the central oscillator recovers the dynamics of a free oscillator. It rather 
leads to a dilatation of all time scales. For instance, the renormalized oscillator frequency 
( vanishes even faster than the inverse damping time, since ( = u / (1 +7 2 ). An analysis 
of the correlator in Eq. (fT2*|) shows that for i < r the particle behaves as a free ballistic 
(though very slow) particle, Cqq (t) oc tj 



T. 



Instead of using criterion A, we may focus on D q (u) = ImCqq (ui)/u. As a spectral 
function, lmC qq \uj) is independent of temperature. Of special interest is the slope of 
D q {uj) near u = 0. Since lim^oo D q {oj) = 0, the condition for the existence for D q {uj) 
displaying a maximum can be written as 

D' q (0) > Criterion B , (19) 

which may be viewed as an indicator of underdamped oscillations or, equivalently, 
coherent transitions. It is often employed for the spin-boson model [231 EHl- For Ohmic 
damping, one finds 



D ^ u > = T7TT r~2 212T 7 7 xT^ > ( 20 ) 



and the critical curve for 7 p is thus given by the relation 7p nt = 7^ (jq/rj 2 — 2) /n 2 , so 
that criterion B is satisfied for 7 P > 7™*. 

Both being based on exact expressions, we observe that criterion B differs 
substantially from criterion A, see Fig. ^ Surprisingly, even for the g-oscillator there is 
a region \/2 < 7^ < 2 where criterion A and criterion B are different. In Fig. ^D q (cu) is 
plotted for different coupling strengths j q , j p . For large couplings j p , we observe that, 
although we are (according to criterion B) in the underdamped region, the maxima of 



1.2 




Figure 2. The spectral function D q {u>) for different coupling strengths j p , j q . On 
the left, the marginal case of the g-oscillator, at (7^, j p ) = (y/2, 0), (- - - -), develops 

a maximum for (j q ,j p ) = (s/2, \/2), ( ). On the right, for a value of 7 9 = 3 

(overdamped regime), when the values of 7 P are ( - - - -), 21 ( ) (marginal 

case), and 40 ( ). 
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symmetric 


only 7g 


only 7 p 


K < 1 


(A) 


always 


lq < 2V 


7p < 2rf x 


D q (u) has peak 


(B) 


7 < 




always 


|Imu;±| < Reu;± 


(C) 


7 < 1 


7 g < y/2n 


7 P < v 7 ^ -1 



Table 1. Condition for coherent dynamics according to three different criteria (left 
column, labeled A, B, C), for three particular cases (upper row). The symmetric limit 
includes the assumption rj = 1. General case is given in main text. 



the curves are not very pronounced. This is related to the fact that the renormalized 
oscillator frequency ( vanishes for large 7 P faster than the inverse damping time r _1 . 

The foregoing discussion suggests the introduction of a third, more restrictive 
criterion for underdamped oscillations, namely )2~Tj . 

C > r- 1 Criterion C. (21) 

According to this condition, the region of underdamped oscillations is convex, i.e. 
by increasing the coupling strength no transition from overdamped to underdamped 
oscillations can be realized. Table I summarizes the different criteria for coherent, 
underdamped dynamics as they are applied to three prototypical cases. 

For high temperatures (fcgT ^> uj ) the integral for the coordinate autocorrelation 
function Eq. (fT2"|) can easily be solved yielding the classical solutions of the Heisenberg 
equations of motion. For k < 1 they read 



exp (^L) • (22) 

For k > 1, Eq. (|22|) describes the corresponding non-oscillating solution. Analytical 
expressions for &+\t) are also available for T = 0. In this case the integral in Eq. (|T^|) 
can be expressed in terms of exponential integrals ^H]- We do not show the result here 
and limit ourselves to note that for long times only the lowest order term, i. e. the term 
linear in u, contributes in the numerator of the r.h.s. of Eq. (|12p. Thus the long time 
behaviour is the same as for the g-oscillator, namely, we have 

2 

C '& ) (*) = -^^ + 0(r s ), for t^oc. (23) 

For the classical solution Eq. (}2"2"j) criterion A [Eq. 1)16)1] is the only natural criterion to 
differentiate between underdamped and overdamped oscillations, since it distinguishes 
solutions with infinitely many zeros from solutions with no zeros at all. For zero 
temperature the number of zeros of Cq^ is always finite and criterion A is less 
informative. 



2 
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The zero temperature mean squares (q 2 ), (p 2 ) can be calculated exactly. They are 
given by 



<s 2 > - m 



4k (1 + 7 9 7 P ) 



^ + 7 p (1-2k 2 ) 



+ — (i + 7p7jin — 2 +o{n ,n ), 

7T CJ 

n, S 1 , K + — 1 , \ 

/(«) = 7=a=F In y= • (24) 

The corresponding expressions for (p 2 ) are obtained by interchanging p <-> g in Eqs. (|24j) 
and (JT3J). For small 7 g ,7 P we have f(n) = 1 — 2k/ it + C(k 2 ) and the position mean 
square becomes 

« = ^-d + ^( ln S-0 +Ofe *'' (25) 

and, correspondingly, 



(P 2 ) = \-^lv + % ^ - 5) + 0(7 9 , Iv) , (26) 

For 7 9 ^ 0, 7 P = we recover the results for the g-oscillator, with the characteristic 
logarithmic dependency of (p 2 ) on the cutoff JH]- I n the general case, the Heisenberg 
uncertainty product diverges as (p 2 )(q 2 ) oc \iaQ p \nfl q . Thus the reduced density matrix 
at equilibrium becomes approximately the identity, its off-diagonal elements being 
essentially zero in both the position and the momentum representation. 

2.3. Other spectral Densities 

For general spectral densities a q and a p are arbitrary positive real numbers. While we 
always have 

lmJ n (uj) = -aga(u)J n (\u\) , (27) 

we find 

{0 , if a n < 2, 

^T>§; if« n = 2, . (28) 

First, we briefly recall the behavior of (q 2 ) and (p 2 ) as functions of a q and 7 g in the 
case of the g-oscillator. For (q 2 ) it may be summarized by the formula jTHj 

d{q 2 ) | < for a q < 2 
d'jq 1 oc — ^ for a q > 2 

A measurement of the position, i.e. a diagonalization of the density matrix in the 
position basis takes place only for a q < 2. The behaviour of (p 2 ) is opposed to that 
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of (g 2 ) in such a way that the product {q 2 )(p 2 ) is always > 1/4, as required by the 
Heisenberg relation. For (p 2 ) we may write 



d{^)_ f oc In (g) for a q < 2 
^ > for a„ > 2 



(30) 



In the forthcoming discussion we focus on (q 2 ). We have a contribution from J q (u), 
which reduces (q 2 ), and another contribution from J p {ui) enhancing (q 2 ) [see Eq. (J55j) ]. 
From the results for the g-oscillator we expect that the J p (w) contribution should always 
dominate for a p < 2 and the J 9 (o;) contribution should be negligible for a q > 2. The 
opposite regime a q < 2 and a p > 2 is more interesting (the case a q = 1 and a p = 3 has 
been studied in Ref. j2I]) • There we have a sum of two terms of order 1 in Q~ l with 
opposite sign. In the following we focus on that case. 

We assume the susceptibility x( u ) to be an algebraic function. This entails 
rational exponents a n = (3 n /m with integer /3 n and m. The integrand in Eq. f!12l) 
has r = max.(/3 q + (3 P , 2m) poles on the m Riemann sheets which split into two classes. 
For aig, a p > 1 there are 2m poles close to the oscillator frequency ±Uq on each sheet. In 
addition there are = max(/3 g + (3 P — 2m, 0) poles of the order of ifl. If either a q or a p 
becomes equal or smaller one, the 2m poles detach from ±uq moving into the complex 
plane. The equilibrium mean square can be written as a sum of the contributions from 
the two types of poles. 

(q 2 ) = C„ + C n • (31) 
For the second term we find (7,, 7 P small) 

Cn = ^^—Q(a p - 1) + 0(rf p , lq ) . (32) 



For a sketch of the derivation of Eq. (}3*2*|) see Appendix A Therefore, the high energy 



poles yield, to first order in 7 p , a finite contribution to the position mean square. This 
contribution is positive for superohmic coupling to the momentum but vanishes for 
Ohmic and subohmic coupling. It has a singularity for a p — 1 which marks the transition 
to the logarithmic dependence on the cutoff frequency, see Eq. (|2"I|). We note that the 
right hand side of Eq. (}32|) does not depend on the oscillator frequency u , i. e. it does 
not depend on the properties of the system itself. The contribution from the poles close 
to ujq, C UQ) can in principle be calculated in a similar way as for arbitrary exponents 
a q , a p ; however, the calculations become increasingly messy. A general treatment is 
also hampered by the wide variety of casuistic behavior, with different regimes defined 
by the conditions a n ^ 1, a n ^ 2 and a q + a p ^ 2, a q + a p ^ 4. 

We focus on the specific case a q < 1 and a p > 2, which is the most relevant, since 
most baths occuring in nature are either Ohmic (Markovian approximation, electron 
gas or superohmic (photon or acoustic phonon baths). In that case, the low energy 
poles are essentially determined by the lower spectral exponent (here, a q ). This reflects 
the general property of dissipative quantum systems that an environment becomes 
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Figure 3. Left: Zero temperature mean square (q 2 ) as a function of 7 9 and 7 P for 
a q = 1 (Ohmic) and a p = 3 (superohmic). Right: Zero temperature mean square 
(<? 2 ) as a function of 7 g for a q = 1 , 7 P = 1 and different spectral exponents a p = 3 

(- - - -), a p = 2.2 ( — • — ) and a p = 1.8 ( ). The full line corresponds to the 

q-oscillator 7 P = 0. The cutoff frequency is Q,/loq — 200 in both figures. The case 
a p = 2 is avoided because of its singular character [see Eq. I|28|) ]. 



increasingly efficient with decreasing spectral exponent a n . To leading order in 7 g , 
we find 



a 



1 

2 



7 9 



- + 3a q (a q - 1) 



+ 0(7, 



(33) 



If we compare C^ with Eq. (|3*2"|) . we observe that the two contributions have opposite 
sign. In particular, the contributions of C uo and Cq cancel each other to first order in 
7 g and 7 P provided that 



1p 



a p — 1 



i + 3n 



9 V^q 



7, 



(34) 



In Fig. 12.31 (g 2 ) is plotted against the coupling strengths 7 P and 7 g for a q = 1 and ct p = 3. 
(g 2 ) is monotonously increasing as a function of 7 P . For 7 p = it is a monotonously 
decreasing function of j q . For the parameters chosen in Figure 12.31 Equation 
holds on the diagonal 7 g = j p . There (q 2 ) remains close to its unperturbed value 1/2 if 

7g,7 P < 1- 

Finally, on the right hand side of Figure EH3 the position mean square (q 2 ) is plotted 
as a function of 7 9 for different spectral exponents a p . For small 7 9 the enhancement 
of (q 2 ) due to the coupling to the momentum is larger for smaller values of the spectral 
exponent a p . This is the expected behaviour. This behavior, however, is inverted when 
7q increases. Then the relative enhancement of (q 2 ) as compared with the g-oscillator 
is bigger for higher spectral exponent a p . A crossing occurs always for some value of 
the coupling constant 7 ? , although some times this can be very large. 
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3. Time evolution 

To study the nonequilibrium properties and in particular the loss of coherence for an 
initially pure state, the calculation of Cqq(t) is not sufficient. Instead, one should 
perform a full nonequilibrium calculation with inclusion of specific initial conditions by 
means of a Laplace transformation. The position-momentum symmetry of the oscillator 
Hamiltonian suggests the use of the reduced Wigner-function W(q,p,t) instead of the 
reduced density matrix. It was shown in Ref. [HI] for the q oscillator that the time 
evolution of W(q,p,t) can only for certain initial conditions be described by an exact 
master equation. Here we consider decoupled initial conditions Wq = WsW^Wg, 
which fall into this class. In this case we are able to derive not only an exact 
Master equation for the reduced Wigner function but a solution in terms of a two- 
fold convolution integral. 

Generally W(q,p,t) can be expressed as 

W(q,p,t) = (8( q -q(t))6(p-p(t)))o , (35) 
where the bracket denotes the average over initial conditions (TTj 

(...)o= / dq'dp'Y[da qk da pk ,da* qk da* k ,(. . .)W s (q',p', {a g , a p , a*, a*}) , (36) 

k,k' 

and q(t),p(t) are the classical solutions of the equations of motion Eq. (JOJ). The Wigner 
function of the thermalized baths is given by 

wW = TT ! exo ( 2a " fcQnfc ^ f37) 

W-7rcoth(/?u; fe /2) V coth(J3u k /2)J ' 1 ' 

the system itself being in an arbitrary pure state characterized by Ws- In this case one 
can express Eq. (|33|) as a twofold convolution integral 

/dk dk' 
™*k P+l k qWs{kUt) + k>(j) _ 
111 Z7l 



exp {-k 2 X pq {t) - k' z X p0 {t) + 2kk'Y pq (t) 
exp (-k l2 X qp (t) - k 2 X q0 (t) - 2kk% p (t)^j , (38) 

where Ws is the Fourier transform of Ws in both arguments. We have introduced the 
auxiliary functions 

<Pn(t) = J dux(uj) (u n + J n {uj)j e lu>t , n = q,p 

Mi) = / duux^y^ , (39) 
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and the temperature-dependent quantities 

1 



X pq (t) 

Xpo(t) 
Y Pg (t) 



2 

1 
2 
1 
2 

Re 



dt(j) q {t -t')e 



dtcj)o(t — t')e 



duJ p (uj) coth(/3a>/2) 

dujJ p {uj) coth(/3a;/2) 
duiJ p {uj) coth(/3a;/2) 
dt(p (t - t')e luJt / dt<f) q {t - t')e 



iujt 



-iujt 



(40) 



The functions X g0 , X^ and y gp are defined accordingly by interchanging p and g. For a 
better understanding of the forthcoming discussion in Sec.|3]we notice that the functions 
defined in Eq. are the time-dependent coefficients of the solutions q(t), p(t) of the 
initial value problem Eq. (JUJ). For instance 



q(t) = Mt)v(0) + </>p(t)p( Q ) + / Mt - t')F p (t')dt' + / <f> p (t - t')F q {t')dt' . 



(41) 



This necessarily requires 0o(O) = 1 and n (O) = 0. Finally, we note that W fulfills a 
Fokker-Planck type of equation |T5] . 



W(q, p, f ) = V [g(q, p, t) + B(t) V] W(q, p, t) , 

where V = (d/dq,d/dp), with g = (g q ,g p ) the phase-space drift term, 
g n (q,p,t) - - ^ G nm {t)m , 



G(f) 



m=g,p 

00^0 + 0p0<? 
0O0q - 9 0O 



^wo - </Wp 

^000 + 0g0p 



+ 



(42) 



(43) 



(44) 



ro 1 v<?rp 

and with (the 2x2 matrix) B(t) the state-independent phase-space diffusion term, 

Bp P (t) = X pq + X ?0 — 2F pp (Xpq + X q0 ) + 2F P9 (F P9 — Y qp ) (45) 
B qq {t) = X qp + X p o — 2F qq [X qp + X p0 ) — 2F q p(Y qp — Yp g ) 

Bpq(t) = Yqp — Ypg — (-Fpp + F qq )(Y qp — Ypg) 

— F pq (X p o + X qp ) — F qp {Xpq + X g ()) ■ (46) 

Bqp{t) is obtained from B pq (t) by exchanging g and p. To derive the coefficients (j4^|) and 
(|45|) we have employed the ansatz Eq. (|%2")l for the Fokker-Planck operator. Acting with 
yet undefined functions f n ,B nm (n,m = p,q) on the r.h.s. of Eq. (JI2J), and comparing 
with its time derivative, yields conditions for the functions f n ,B nm . An exact Fokker- 
Planck equation in the form of Eq. (1421) for the g-oscillator has been first derived in 
Ref . jTT] and later by a different method in Refs. [T3J EH] . 
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3.1. Purity 



A convenient quantity to measure the degree of global decoherence is the purity 
V(t), defined as the average of the density matrix itself [T7] . 

P{t) = {p{t)) = trp 2 (t), (47) 

which is basis independent. Of special interest is the equilibrium value Vp = lim^oo V(t) 
which measures the efficiency of the environment in destroying quantum coherence. Here 
we implicitly assume ergodic behavior which, as we shall see, applies in the presence 
of Ohmic baths. For a harmonic oscillator in thermal equilibrium, the reduced Wigner 
function is [TH] 

Wflb'P) = 7, , l1/2 ex P ~ tAA , (48) 



which leads to 

Vp = or/ 2^ h i 2^ 11/2 • ( 49 ) 

4. Coherence decay for Ohmic damping 

Although the equilibrium decoherence (as measured by the product (q 2 ) p(p 2 ) p) is 
enhanced by the additional noise term, one may wonder whether for low temperatures 
the decoherence time becomes larger than for the damped g-oscillator. To answer this 
question exhaustively one would have to calculate the time evolution of the purity for 
an arbitrarily pure initial condition. We consider two cases. First, we choose a coherent 
(Gaussian) state as initial state. This case should present the greatest robustness against 
decoherence [T71 . Second, we choose a superposition of two Gaussian wave packets. 
Then a new aspect of decoherence comes into play, namely, the fast vanishing of the 
relative coherence between the two Gaussian wave packets. We distinguish the two 
manifestations of decoherence by introducing a new quantity, which we call relative 
purity. 

4-1. Decoherence for a Gaussian wave packet as initial state 

We consider here the case where the system starts in a coherent state at t — 0. The 
Wigner function for a Gaussian wave packet is 



W s (q,p) = -exp 

7T 



2 1 / \2 



-v (q -qo) - - (p - Po) 
v 



(50) 



where p = y^Im (a) and q = (a) are defined in terms of the complex 

eigenvalue a (amplitude) of the coherent state (recall r\ = ^Juo q /uo p ). Since Ws(q,p, 0) is 
Gaussian, the convolution integral of Eq. (|38|) can be performed. The integrals become 
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particularly simple when both baths are equal (7 P = 7 g = 7 and 77 = 1), corresponding 
to the completely symmetric case. Then we have <ft p = <p q = <\>\ and, by the same token, 
X pq = X qp = Xi in Eq. (jHEJ) and also X p0 = X q0 = X in Eq. (jjQ)) . The crossed terms 
in Eq. ([38)1 vanish and the Wigner function becomes a product of two Gaussians for all 
times 

W(q,p,t) 1 ' 



4vr^ + 0f + 4X + 4Xx 



exp 



(p - Po0o 



+ {q 



+ 



o : + # + 4X 1 + 4Xo 
For the purity one obtains the simple expression 



P(f) 



1 



$ + # + 4X„ + 4Xi 



(51) 



(52) 



The remaining task is to calculate the quantities 0o ; X , Xx from their definition in 
Eqs. and (f4f)|) . Using the naive form (Q — > 00) of the spectral function of an Ohmic 
bath as given in Sec.El we would find 0o(O) = 1/ (l+7 P 7 g ) + 1 1 leading to inconsistencies 
[see the discussion after Eq. (jlHj) ]. The reason for this lies in an initial slippage caused 
by the somewhat unphysical character of the decoupled initial condition jHSl EH EEl IHHj ■ 
Technically it stems from the fact that for J(u>) oc u the integrals in Eqs. (|3~9j) do not 
converge at t = This problem is overcome by regularizing the Ohmic spectral 
functions, i.e. by reintroducing a finite cutoff Q. The explicit calculation with a Drude 
regularized spectral function is sketched in |Appendix C| Here we state only the main 
result: At T = the purity decays on two time scales, given by Q^ 1 and r in (fIH| . 



V(t) 



V* 



1 



o^t 5 i+r 



e 



■nt 



■ cos(At)e 



-i/T 



for < t < fi- 1 
for < t -> 00 



(53) 



with P^" 1 = 2(q 2 )p and the oscillator frequency A = co>o/(l + 7 2 )- Coherence is reduced 
immediately after the start of the coupling. Although afterwards it decreases more 
slowly, on a time scale r, for larger couplings the curves of V(t) for different values of 7 
never cross, i.e., V(t) is a monotonously decreasing function of 7 for all t. 

An initial slip similar to that discussed in this section also occurs for the g-oscillator 
[11 . However in that case its effect on purity is much less severe. Specifically, out of 
the set of functions (fi q (t), 4> p (t) and (f>o(t) defined in Eq. ()39|) . only (f) q (t) is affected by 
the initial slip. After a time it becomes <f> q (t) ~ 7 P + 1 JT]. However, we note that 
4> q appears in the temperature-dependent terms X qp , X pq , Y qp , Y pq of Eq. (j4Ti|) only in 
combination with the spectral density of the momentum coupling J p (uj), which vanishes 
for the g-oscillator by definition. We thus reach the important conclusion that the 
purity evolution of the g-oscillator is insensitive to the initial slip stemming from the 
use of decoupled initial conditions: At t ~ Q^ 1 the purity is still approximately unity, 
decreasing afterwards at a rate oc j q . 

Isar and coworkers [10] studied in detail the decay of purity for the g-oscillator. In 
particular they found constraints that must be satisfied by the bath if the purity is to 
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remain constant and close to unity during the whole time evolution of the oscillator. 
Such a high purity is not possible in the present case due to the initial slip [see Eq. (153] . 

Another possible preparation of the initial state is in a constrained equilibrium. The 
expectation values (q(0)) = qo and (p(0)} = p are held fixed with the bath equilibrated 
around those values. Then the Wigner function is ^1] 

w {q , P ,t) = }-—^(J<-W)?_<>-W>)?\ (54) 

K " ' V 2{, 2 ) 3 2(p>)f, J 

while the purity is given by its equilibrium value Eq. (|49|). A more detailed discussion 
of the purity decay of the g-oscillator, together with the general formulae, is given in 
Appendix C| 



4-2. Decoherence of two Gaussian wave packets 



Formula (]3q) also allows us to investigate more complicated initial conditions such as as, 
for example, the superposition of two Gaussian wave packets. This case has been studied 
for a single bath by Caldeira and Leggett [H]. It is an interesting case study because it 
displays two different aspects of decoherence. On the one hand, there is the decoherence 
which either wave packet would experience alone. This part is essentially described by 
V(t), which we will call Gaussian purity [see Eqs. (152]) and (|53)l ]. On the other hand 
there is the decoherence due to the spatial separation of the two packets. This second 
contribution is expected to become increasingly important when the distance a between 
the two packets becomes large. As initial wave function we choose 

, . . 1 / (x + a/2) 2 (x-a/2) 2 \ 

^ = 4/7^ e + e > ( 55 ) 

c\l Iko 1 \ J 
which translates into an initial Wigner function 

W s (q,p) = \e-^^ J2 + e ~ lkaP ) ■ ( 56 ) 

Here c = y/2 [1 + exp(— a 2 /8cr 2 )] 1 ^ 2 is a normalization constant. We assume the most 
symmetric case and set a = l/y/2. Plugging (|5fi]l into (jHHj) we find that the purity can 
be expressed in terms of the Gaussian purity as 



V{t) 

The function <fr(t) is given by 



cosh 



p(t) = ^< 1 + — cosh- ( aV8) > (57) 



where 0o an d 4>i are defined in (J39"]) and computed in ()B.2|) for the symmetric case. <p(t) 
evolves from 0(0) = 1 to lim^oo <p(t) = 0. We define the relative purity as the ratio 

Prel(t) = P{t)/V{t) (59) 
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As expected, P Te i(t) — > 1 as a — > 0, and P re i(t) — > 1/2 as a - >oo. 

Interestingly, the structure of (|57jl is such that, as time passes and <p(t)V(t) evolves 
from 1 to 0, the ratio P re i(i) starts at unity, as corresponds to a pure state, then decreases 
and finally, at long times, goes back to unity. When a is large P Te \(t) decays rapidly to 
1/2, on a timescale ~ l/4a 2 7. There it stays for a time which increases with distance 
as ~ 7 _1 lna. Afterwards it returns to one. The ratio 1/2 can be rightly interpreted 
as resulting from the incoherent mixture of the two wave packets. Thus it comes as a 
relative surprise that P re i(i) becomes unity again at long times, as if coherence among 
the two wave packets were eventually recovered. The physical explanation lies in the 
ergodic character of the long time evolution, with both wave packets evolving towards 
the equilibrium configuration [see Eq. and the subsequent discussion]. 



5. Single bath with linear coupling to both position and momentum 

It is instructive to compare Eq. © with the Hamiltonian of a harmonic oscillator 
interacting with a single heat bath in the most general form of linear coupling, 



a k H q H P 



(60) 



with complex parameters \k, l^k- This model is different from that discussed in the 
previous sections in that here time reversal invariance is broken. By this we mean 
the following. In the models described by Eq. (JHJ), the bath modes might describe 
for example a magnetic field coupled to the momentum of a charged particle, which 
clearly would break time reversal invariance. However in that case, such a symmetry 
breaking is somewhat fictitious since, due to the linear nature of the coupling and to 
the modelling of the bath as a set of harmonic oscillators, one can always find a unitary 
transformation which restores time reversal invariance, i.e. which renders all parameters 
in the Hamiltonian real quantities. It is easy to see that, for general (complex) and 
fik, such a unitary transformation cannot be found for Eq. (}60j) . 

By an analysis similar to that of Section 121 one finds general expressions for the 
symmetrized equilibrium correlation functions: 

C ( +\t) =-/ cMxM| 2 cos(^)coth(fy3u;/2) 



7T 



+ lmJ p (u>) 



ImJ 9 (a;J 



u v + Re JJu) 



u + ReJ-(Lu)) + (Re«7+(w) 



uj p + Re J p (u>) 



Im 



2uJ_(uu) 



Jl(u) 



(61) 



The same expression applies for Cpp\t) with the indexes p, q interchanged everywhere. 
As before, we define a generalized "susceptibility" for the system which now reads 



LUq - J q (u) 



^ P — J P (u) — u 2 — J+(v) — J^i 1 ^) + 2cjJ_(cg>) 



(62) 
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Eqs. and (|62|) involves four different spectral functions. J g (u) and J p {u) are defined 
as in Eqs. (@J) and (jlOj) . Now we introduce the spectral functions 

= 2^|Ajfc||jt/jfc| cos9 k S(uj - u k ) 

J-(w) = i J^(A fe /ifc - \* k fx k )5(uj - u k ) 

= 2^ |Ajfc||jU fc | sia9 k 8(u; - u k ) , (63) 

which reflect the mixing of the bath modes. The "hat" symbol denotes the 
transformation 

f(co) = uV r ~^- 2 dw' + £f(W) , (64) 
J uj /z -uj 2 2 

In Eq. (164)1 the real part is an antisymmetric function and the imaginary part is 
symmetric. This is exactly reverse to the "tilde" (Riemann) transform, defined in 
Eq. (HHJ). 

Before we calculate Eq. (}6*Tj) for a specific case, we analyze some generic features. 
First we notice that, by setting J-(ui) and J+(u)) to zero, we recover the autocorrelation 
function for two independent baths [see Eq. ()12|)]. On the other hand, the two spectral 
densities J+(lj) and J- (a?) are not independent from J p (u) and J q {oS)\ rather, they 
satisfy 

J 2 + {uj) + Jl(u) = J p {uj)J q {uj) . (65) 

Actually this relation has already been used to simplify the integrand of Eq. (|61j). We 
note that it holds for the spectral densities themselves but in general not for their 
Riemann transforms J+{oj) and J-(uj). Apart from the above constraint one can, at 
least in principle, freely choose three of the four spectral functions. If, in addition, 
we make the physically reasonable assumption that the "mixing angle" in Eq. (|63)) is 
frequency independent (8 k = 9), condition fixes J+{oj) and J-{oj) completely. 

J+(uj) = \J q {uS) Jp{w)\ COs9 

J_[u) = [J q (u)J p (uj)} 1/2 sm9 . (66) 

Therefore, in the important case where both J q (uS) and J p (uj) have the same spectral 
exponent (a q — a p — a), J + and J_ also obey the same power law. Using Eq. ()63)) we 
observe that the term (Im J g Im J p ) implicit in Eq. ()62j) drops out. We recall that this 
interference term has been responsible for the non trivial "phase-space" diagram Fig. 
UJ What is more, only Re J + appears in Eqs. (|6H]l and (|52|) . That means, according to 
Eq. (J2EJ), that for a q + a p < 4 the symmetric spectral function J + does not contribute 
at all. The relation between the "hat" and the "tilde" transformation is 

f{u) = J/M , (67) 

from which all properties of / can easily be deduced. For example from Eq. (|67)1 and 
Eqs. (J2HJ) and (j23) it follows that ReJ_(cu) = only for a q + a p < 2. 
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It is illustrative to look at the "generalized" susceptibility Eq. (J52j) for the case 
a q + a p < 2, where all real parts of the spectral functions vanish. We obtain 



X 



Uq — {J 2 — uu q J p (iju) — u p J q {uj) + 2uJ-(u) 



uj — on — sgn [u! ) 



2\uj\ 



J g (H)+^J p (H)--L-Lj_(H) 



Mr, 



(6* 



This is exactly the suceptibility of a g-oscillator coupled to a single bath with an additive 
spectral function 



J(uj) = J q {uj) + rj 2 J p (uj) + 



2uJ-(u) 



0J„ 



(69) 



We conclude that, when only one bath is involved, the particular structure of the 
coupling of the bath to the position or momentum variable can always be modelled 
by a g-oscillator with an appropriately chosen spectral function. In this context the 
combination uJ_{&) can be considered as stemming from an effective additional noise 
source that couples to position. Note however, that this is not the same as saying that a 
single-bath coupling to q and p can always be unitarily transformed to a coupling only 
to q (^-oscillator), since this is not true in general. For one thing, the noise properties, 
which are not uniquely given by x( w ); would be different. 

We also wish to emphasize that a double-bath dissipative oscillator cannot be 
described in terms of an oscillator coupled to an effective single bath. In particular, 
it can never be modelled by a susceptibility like (|69jl. This is the reason why the physics 
explored in Section El is so different from that which could be found in any possible 
single-bath scenario. 

We do not repeat here the analysis of Section El for arbitrary spectral functions but 



focus instead on the case of Ohmic coupling (a q = a p = 1 
the results obtained for two independent baths. As Re J n (a;) 
Section l2~3j . we obtain 



in order to compare with 
= ReJ+(co>) = vanish [see 



7T 



e?u;|x(a;)| cos(u;£) coth.{hj3uj / 2) 



Im 



uj q J q (uj) + J p (uj) 



id 



ReJ(u) 



2ui p J_ 



uj — Re J-(u>) 



.(70) 



The form of the susceptibility varies depending on whether the mixing angle is a multiple 
of 7r or not. We discuss the two cases separately. 



5.1. Mixing angle 6 = 0, it 

For 9 = or 9 = n one can bring, by an appropriate redefinition of a^, a k , the interaction 
part of the Hamiltonian Eq. (J60|) into the form 

E I oc ^q k {W\p+ |A fc |g) , (71) 

so that the bath couples to the main oscillator either through the position of their 
oscillators or through their momentum, but not through both. From Eq. (|6l)|) it can 
be seen that the antisymmetric spectral function vanishes identically. Apart from the 
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Figure 4. Regions of miderdamped oscillations for a single bath coupling to position 
and momentum (dark region) and for two independent baths coupling to position and 
momentum (dark and bright region). r\ = 1. 



susceptibility, the integrand of Eq. (|70|) becomes equal to that of Eq. (pH^j) . Thus all 
results of Sec. 12.21 m particular Eqs. (|22|) and (|24|) . carry formally over. The only but 
crucial difference lies in the poles of the generalized susceptibility x{ u )- The set of zeros 
can be written as 

u± = too (-zk/2 ± Vl - k 2 ) , (72) 

with k defined as ^ 

K = I 7 ^ + 2^ • (73) 
According to criterion A, Eq. (fTtjj) . the regions of overdamped and underdamped 
oscillations in the (jg, j p ) plane are separated by a straight line with negative slope given 
as before by the condition k = 1 [see Fig.|3]. Therefore, a crossover from overdamped to 
underdamped oscillations by increasing one of the coupling constants is now impossible. 
We also note that, in the symmetric case (7^ = 7 P = 7), the relaxation time is the 
same as for the g-oscillator [r _1 = uj j/2\. Increasing 7 g or 7 p leads inevitably to an 
enhancement of dissipation. This is in blatant contrast to the results obtained previously 
for the double-bath model and one of the most remarkable results of this work. It best 
illustrates the importance of the specific structure of the bath and justifies a posteriori 
the detailed study of a model with two independent baths. 

5.2. Mixing angle 9 ^ 0, 7r 

We briefly consider the case of non-vanishing mixing angle. The effects should be largest 
for 9 = 7r/2, corresponding to an interaction Hamiltonian 

Hi oc } j (\Hk\PkP+ \h\qkq) ■ (74) 

This form of coupling is called amplitude coupling in |43J. Now J-{oj) is not longer zero 
but we have 

luJ-(uj) = A /7 9 7 P o; 2 



2 / Q 2 
In ( H ) + isgn (u) 



71 



(75) 



Dissipative quantum oscillator with two competing heat baths 



22 



which is exactly the spectral function of a superohmic bath with exponent a = 2. The 
logarithm in Re J_ inhibits further analytic treatment of Eq. ()7U|) even in the high 
temperature limit; thus we limit ourselves to a qualitative analysis. For simplicity, we 
omit the terms containing logarithms in the generalized suceptibility Eq. ()62p. The 
contributions to Cqq stemming from those terms should decay on a time scale ~ O -1 . 
The remaining terms define a suceptibility which is identical to the suceptibility of a 
g-oscillator coupled linearly to one bath with spectral density 



It is well known that, in the case of a polynomial spectral function, the term with 
the lowest exponent dominates [see also the discussion in Sec. 12. H] . Therefore we may 
conclude that the crossover diagram Figure 0] remains unchanged by a mixing angle 



6. Conclusions 

We have discussed the behavior of a quantum Brownian particle in a harmonic potential 
subject to two independent noise sources, one of which couples to its position and the 
other one to its momentum. In the symmetric case where both baths are Ohmic and 
their coupling strength is the same, we find underdamped oscillations of the central 
oscillator for all coupling strengths. This indicates that the two baths partially cancel 
each other. The effect is due to the mutually conjugate character of position and 
momentum. It was first noted in Ref. [2E] for the (cylindrically) symmetric spin-boson 
model with spin s — |, i.e. in the deep quantum regime. "Quantum frustration" of 
the spin can pictorially be described as the result of having two observers attempting 
to measure simultaneously, with equal efficiency, two non-commuting components of 
the spin. Because of the uncertainty principle, both of them fail to measure anything. 
Here we have investigated the analogous effect for a quantum oscillator coupled to two 
independent baths. We have found a moderate form of cancellation which we have 
labelled "quasiclassical frustration" because our dissipative oscillator may describe a 
large spin impurity coupled to the magnon bath in a ferromagnetic medium. The most 
notable features of quasiclassical frustration become manifest in the strong coupling 
limit, where underdamped dynamics survives and all time scales diverge. It remains to 
be investigated whether the occurrence of frustration in the classical regime is a general 
property or an artifact of the harmonic oscillator. 

We have compared the double-bath model with the case where a single bath 
couples linearly to the position and momentum of the oscillator. In the latter case the 
"phase-space" diagram is simpler in that transitions from overdamped to underdamped 
oscillations can never occur. Comparison of the two models indicates that bath 
correlations can qualitatively change the behaviour of a dissipative system. 

A point of caution is needed in the interpretation of our results. We have definitely 
ruled out the at first sight enticing but at closer inspection unphysical conjecture that 




(76) 
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the effects of the two observers cancel each other completely and the particle is not 
affected by the environment. Indeed we have seen in the specific example of a decoupled 
initial state that, in destroying quantum coherence, two baths are always more efficient 
than one. This is true for both the overdamped and the underdamped regime. This 
demonstrates that, at least for the harmonic oscillator, underdamped dynamics is by no 
means a reliable signature of high global coherence, which here we identify with purity 
[see Sec. 0J. Our work provides further evidence that decoherence and dissipation are 
not necessarily correlated. By dissipation we mean here the net transfer of energy from 
the central oscillator to the thermal baths. It is characterized by the classical equation 
of motion or, equivalently, by the properties of the spectral function. We have seen 
that, depending on the situation, an increase of dissipation can be accompanied by 
either a reduction or an increase of decoherence. Vice versa, a source of decoherence 
may or may not lead to dissipation. This is the case in the so called pure dephasing 
(not considered here), where the interaction part of the Hamiltonian commutes with the 
system Hamiltonian and no dissipation occurs at all, see i. e. |43j . 

The fact that we have focused on a quantum oscillator coupled linearly to 
oscillator baths has allowed us to investigate analytically the equilibrium and dynamical 
properties. The question of frustration is also raised in other dissipative quantum 
systems, such as a small spin coupled to a boson [2B] or a fermion jUj bath, which 
are not amenable to an exact analytical treatment. The extension of the present study 
to less tractable physical scenarios provides a theoretical challenge. 
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Appendix A. Calculation of Eq. (1321) 

We calculate the high frequency contribution to the position mean square in Eq. (}3*T]) . 
This contribution vanishes for a q + a p < 2 but yields a finite contribution otherwise. 
We assume again a q = f3 q /m and a p = P p /m, (j3 q ,/3 p ,m e Z) to be rational numbers. 
Then this contribution can be written in a determinant form 



Here we use the Vandermonde determinant A^(A 2 ) = rii<j(A? ~ ^|)- We recall that 
r is the total number of poles, while < r is the number of high energy poles. The 
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function f(x) is given by 

40(0; _ 2) 

f(x) = p p + 3m - 1, r) + lqlp v p J 2 $(x, 5m + f3 q - 1, r) 

+ 7,7 p $(x, 2/3 p + P q + m - 1, r) , (A.2) 

where the function <&(x, s, r) is a rather complicated function involving a special function 
called Lerch transcendent [42]. It can, for s odd, be expressed in terms of a logarithm 
plus a polynomial 
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(A.3) 



In order to use this formula we assume m and f3 n odd. Consequently we formally exclude 
a n even. However, that case can be included by considering the limit m —>■ 00 with j3 n 
= 2nm — 1. The argument of f(x) is the square of one of the r roots of the characteristic 
polynomial 
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The asymptotic behavior of Cq for j q , 7 P <C 1 can be derived by expanding $ in Eq. (|A 



For tt p > 1 we find to lowest order 
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7T 2 (« P — 2) 2 a q + 1 2a p + a g — 3 
Keeping only the first term for small 7 ? and small 7 P yields Eq. 

Appendix B. Derivation of Eq. (I53R 

We regularize the Ohmic spectral functions by a Drude cutoff 

27 n VL 2 n U ~ t N 7n^^r 
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Assume Q q = Q p = Q, 7g = 7 P = 7. For 0o and 0i we obtain 
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where no ambiguity is left [A = u; (l+7 2 ) 1 and, we recall, r 1 = 7A]. In a time of order 
Q^ 1 after the connection (with Q — > 00), both functions have dropped on average to a 
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value (1 + 7 2 ) -1 . This is the value one would have obtained by using directly J(u) oc u>. 
After this short initial period, 0o an d <f>\ decay much more slowly, on a time scale r. In 
the expression (JH2|) for the purity decay after the initial slip, we may safely neglect the 
first two terms in the denominator. The functions <po, however still govern through 
Xo, Xi the time evolution of the purity V(t) [see Eq. pOjl]. We notice that both X 
and Xi are zero at t — 0. However in the initial time interval of order I/O, they both 
increase rapidly After the initial slip, they settle to a value ~ (<? 2 )/37 2 / (1 + 7 2 ) and 
increase much more slowly afterwards. This can be seen if we write for instance Xo(t) 
as follows 

X (t) = l [ du J (u) colli {(3u/ 2) IxHI 2 ^ 2 |l - go(uj,t)e- iuJt \ 2 , (B.3) 



where g (u,t) is the principal value integral 



Qofat) = J-^V [ ^p^e^dJ . (B.4) 



27ri x(cu)a; 

Eq. (|B.4J) can be evaluated straightforwardly using the residue theorem. However for 
our purposes the most important point is that g (u,0) = 1 and that go(u, t) behaves in 
time essentially in the same way as (j)o(t) in Eq. (|B.2j) . Thus we may write 

g (cu, t) = — - L cos (At) -i(u + J(u)) sin (At)] e~* /r 

+ [M") cos( 7 2 fit) + B a (u) sin( 7 2 fit)l e~ m . (B.5) 

1 + 7 2 L J 

Here An, Bq are some algebraic functions of uj whose exact form is not important for 

the present discussion, since after the initial slip of time the second line of Eq. ()B.5j) 

can be neglected. Proceeding in the same way with X\(t) we arrive at 

(1 + 7 2 )V l + 7 2 

where c£\t) is given in Eq. p2|l . The temperature dependence is encoded in (q 2 )/3 
and Cqq\t) . Of course the r.h.s. of Eq. (jB.6|) is not zero any more at t — due to 
the neglected initial transient. At zero temperature the autocorrelation function decays 
algebraically oc — u^t' 2 . Then we reproduce Eq. 

Appendix C. Purity decay of the q— oscillator and general formulae 

Here we state the general formulae for the purity decay of a coherent (Gaussian) initial 
state. To calculate the purity decay for two Ohmic baths, both with decoupled initial 
conditions but with different coupling strengths, we proceed as in Sec. 14.11 Performing 
a twofold Fourier transform of Eq. (jHUj) and plugging the result into Eq. (J3%|) yields for 
the reduced Wigner function 



2X {t) + 2X 1 {t) = (q 2 ) p 1 + - ■ 9X9 + - m V cos (At) e"^ , (B.6) 
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+ 11 + 4jc + 4X g0 
77 



A p = ^ + $77 + AX qp + 4X p0 
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2 2r/ 

For the purity we have with the above definitions 



V{t) = {A q A P -AA\ p ) 
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(C.2) 



(C.3) 



As expected, this expression reduces to Eq. (J52J) for J q = J p and r] — 1. For J p = we 
have the case of the g-oscillator. From the definitions (jjUj) we get X p9 = X p0 = Y pq = 
and obtain a closed expression for the purity decay of the g-oscillator, 
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The explicit form of the functions (f>o(t), <f) n (t) with Drude regularized Ohmic spectral 
function can be found in Ref. jllj . Here we give only their form after taking the limit 
Q — > oo, i. e. for t ^> Q _1 . 



cos(Ct) - — sin(Ct) 
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As in the main text ( is the renormalized oscillator frequency. With these functions we 
can calculate X q0 (t), X qp (t) and Y qp (t). We obtain, 



X, 



Y, 



'IP 



>(*) = I 1 + ^Wl - *>(*)C„(f) + \<flLW)p - ^-C£\t) 

(*) = ^ [0o(t)(g 2 ), - <?«(*)] + ^ [( P 2 ),Mt) - c£\t)] (C.6) 



We notice that 



Y qp (0) = X q0 (0) = X qp (0) = P (O) = . (C.7) 

Therefore we find from Eq. (jC.4|) . V(0) = 1, i.e., the purity is unaffected by the initial 
slip. 
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